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A  THERMAL  HYDRAULIC  MODEL  OF  MELT-LUBRICATION  IN  RAILGUN 

ARMATURES 

R.  E.  Kothmann 

Kothmann  Consulting,  Pittsburgh,  Pennsylvania 

F.  Stefani 

Institute  for  Advanced  Technology,  Austin,  Texas 

Abstract — This  paper  describes  the  first  step  towards  a  model  of  the  liquid  film  at  the 
rail/armature  interface  in  solid  armature  railguns.  The  model  considers  high-speed 
Couette  flow  with  purely  viscous  heating  and  does  not  include  MHD  body  forces  or  Joule 
heating.  The  focus  of  the  model  is  on  coupled  fluid  dynamics  and  multi-phase  heat 
transfer.  The  formulation  is  similar  to  the  analysis  of  melt  lubrication  in  rotating 
projectile  bands.  However,  our  first  principles  model  allows  the  possibility  of 
solidification  while  the  armature  is  passing,  a  feature  that  has  been  missing  from  previous 
analyses.  The  model  is  moderately  successful  at  reproducing  results  of  experiments 
which  measured  high-speed  mechanical  wear  of  7075  aluminum  sliding  against  ETP 
copper  for  face  pressures  ranging  from  6  to  22  ksi.  Discrepancies  between  calculated  and 
experimental  results  are  attributed  to  uncertainties  in  modeling  the  complex  phase  change 
behavior  of  aluminum  alloy  7075  and  uncertain  conditions  at  the  rail  interface. 


INTRODUCTION 

The  normal  operation  of  a  solid  armature  railgun  necessarily  involves  melting  and  loss  of 
aluminum  from  the  armature.  This  is  evident  from  the  coating  of  solidified  aluminum  that  is  left 
on  rails  practically  from  the  starting  position  of  the  armature.  Excessive  or  uneven  loss  of 
aluminum  from  the  armature  is  not  desirable  because  it  can  lead  to  arcing  contact.  However, 
there  is  reason  to  believe  that  some  loss  of  aluminum,  sufficient  to  maintain  a  stable,  liquid  film 
interface,  has  several  benefits:  it  provides  an  uniform,  low-voltage  electrical  contact  between  rail 
and  armature,  it  reduces  the  frictional  drag  force  between  armature  and  rail,  and  there  is  evidence 
to  suggest  that  it  suppresses  the  onset  of  gouging  by  several  hundred  meters  per  second  [1], 

There  are  several  processes  by  which  aluminum  is  removed  from  the  armature.  Most  armature 
wear  models  to  date  have  focused  on  “current  melt  wave”  behavior,  whereby  intense  skin-effect 
heating  leads  to  localized  melting  and  loss  of  aluminum  at  the  perimeter  of  the  rail/armature 
interface  [2-7].  When  the  armature  is  moving  faster  than  700  to  1000  m/s,  however,  the  root 
cause  of  wear  is  not  obvious,  since  at  this  point  mechanical  (viscous  heat  generation)  and 
electrical  heating  are  both  potentially  significant. 

The  importance  of  mechanical  heating  or  “melt  lubrication”  as  a  wear  mechanism  in  railguns 
was  demonstrated  in  a  series  of  tests  conducted  at  the  Institute  for  Advanced  Technology  in  1997 
[8].  The  tests  measured  the  wear  rate  of  7075  aluminum  alloy  in  high-speed  sliding  contact  with 
ETP  copper.  In  these  tests,  small  sliders  were  accelerated  from  rest  to  2200  m/s  while 
maintaining  a  steady  contact  pressure  throughout  most  of  the  launch.  Pressures  ranged  from  6.5 


ksi  to  21.8  ksi.  Even  though  the  tests  were  conducted  using  a  railgun,  they  measured  purely 
mechanical  (frictional  and  viscous)  components  of  wear  by  electrically  isolating  the  sliders  from 
the  current  path.  The  tests  found  that,  at  pressures  typical  of  armature/rail  contact,  wear  due  to 
mechanical  heating  was  comparable  to  the  wear  rate  measured  in  armatures. 

The  results  were  originally  presented  as  a  set  of  five  curves  showing  measured  depth  of  wear  vs. 
distance  traveled.  The  same  data  are  shown  in  Figure  1,  plotted  in  a  form  that  is  more  easily 
compared  to  the  output  of  the  model  we  propose.  The  figure  plots  the  velocity  at  which  the 
slider  is  wearing  (Melt  Front  Velocity)  against  the  speed  of  the  slider  in  the  experiments.  The 
data  do  not  extend  below  800  m/s  because,  until  the  sliders  enter  into  a  state  of  melt  lubrication, 
there  is  no  measurable  wear. 


Melt  Velocity  vs.  Slider  Velocity 
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Fig.  1.  Curves  obtained  from  wear  tests  of  7075  aluminum  in  sliding  contact  with  ETP  copper  [8]. 
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DERIVATION  OF  MODEL 


Derivation  of  Equations 


Cu 

constant  in  solution  for  U(Y) 

u,  v,  w 

c,  cr,  cs 

specific  heat  of  melt,  rail,  slider 

u„  u. 

Ec 

Eckert  number 

1*  2 

V 

F, 

normal  force  on  the  slider 

0 

v 

Hm 

latent  heat  of  melting 

ms 

X 

h 

film  thickness 

X 

k.  K  K 

thermal  conductivity  of  melt,  rail, 
slider 

Y 

L, 

smaller  side  of  the  slider 

y 

slider  dimensions  in  x  and  z 
directions 

z 

m. 

mass  melt  rate  per  unit  area 

Z 

n 

index 

P 

8 

Pr 

Prandtl  number 

P 

pressure 

c 

tt>  t. 

temperature  of  melt,  rail,  slider 

0 

incipient  melting  temperature  of 
slider 

k,  k;,  ks 

liquidus  temperature  of  slider 

X 

melting  temperature 

p 

V 

pseudo  melt  temperature 

P>  Pt>  Ps 

x'o 

initial  temperature  of  rail 

a 

u,v,w 

dimensionless  velocity  in  x,  y,  z 
directions 

X 

Table. 

velocity  in  x,  y,  z  directions 
components  of  u 

melt  film  velocity  at  slider  interface 
melting  rate  into  the  slider 
position  of  re-solidification/melt  interface 
direction  of  slider  motion 
dimensionless  coordinate 

direction  X  to  melt  plane 
direction  in  melt  plane  1  to  x 

dimensionless  coordinate 

ratio  of  longer  to  shorter  side  of  slider 

average  pressure  factor 

fraction  of  melt  solidified  under  the  slider 

dimensionless  temperature  in  melt  film 

thermal  diffusivity  in  melt,  rail,  slider 

solidification  rate  parameter 
absolute  viscosity  of  melt  film 
density  of  melt,  rail,  slider 
load  per  unit  area 
time 


Our  model  is  based  on  a  model  of  melt  lubrication  by  Stiffler  [9],  as  shown  in  Figure  2a.  The 
model  assumes  a  block  of  length  Lx  and  width  Lz  moving  at  a  constant  x-directed  velocity  u0  A 
steady  force  normal  to  the  rail  Fy  =  Aa  =  ctLxLz  is  applied  to  the  block.  Viscous  heating  causes 

t  e  s  ider  to  melt  at  a  constant  melt  velocity  vms.  The  film  has  a  constant  and  uniform  thickness 
h. 


Figure  2b  illustrates  the  model  proposed  for  the  railgun  armature.  It  is  based  on  the  same 

assumption  made  by  Stiffler  with  the  exception  that  the  film  material  is  allowed  to  solidify.  This 

is  reflected  in  the  italicized  text  in  assumption  (4)  below. 

1 .  The  fluid  is  laminar  and  incompressible. 

2.  As  per  standard  lubrication  theory,  the  pressure  is  constant  through  the  thickness  of  the  film 
and  average  values  of  density  and  viscosity  are  employed. 

3.  Film  thickness  and  melt  penetration  into  the  solid  are  small  compared  to  the  length  and  width 
of  the  slider. 

4.  Only  one  surface  is  melting.  Solidification  can  occur  on  the  stationary  surface  during  transit 
of  the  slider. 

5.  Quasi  steady-state  conditions  are  reached  in  the  film.  That  is,  we  are  solving  for  a  set  of 
average  conditions  over  the  length  and  width  of  the  contact  face. 
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■Stdtaf,  No  feeling  F  =  Atf  = 


P  ”  ^>>1  tor  problem  of  correlating  weariest 


Fig.  2.  (a)  Schematic  of  melt  lubrication  model  by  Stiffler;  (b)  our  model  is  based  on  Stiffler’s  but  allows 
solidification  under  the  slider. 

Assumption  (3)  is  justified  by  the  experimental  finding  that  the  thickness  of  the  film  is  on  the 
order  of  20  microns  or  less  [8] .  The  film  thickness  of  20  microns  or  less  is  appreciably  smaller 
than  the  smallest  dimension  of  the  slider  (1.4  mm)  and  the  thermal  diffusion  depth  into  the  slider, 
which  is  on  the  order  of  50  microns.  Viscosity  is  approximately  constant  because  there  is  very 
little  variation  in  temperature  across  the  film,  and  furthermore,  the  change  in  viscosity  of  liquid 
aluminum  with  temperature  is  small. 

The  assumption  that  only  one  surface  is  melting  is  consistent  with  evidence  that  rails  are  not 
melting.  On  the  other  hand,  the  assumption  that  the  film  can  solidify  marks  the  principal 
difference  in  bur  approach  and  is  a  key  feature  of  our  model.  Additionally,  Stiffler  neglects  the 
effect  of  mass  injection  on  anything  but  the  squeeze  film  velocity;  we  include  the  effect  of 
injection  on  the  velocity  profile  in  a  way  that  affects  dissipation  in  the  film  and  heat  transfer  to 
the  boundaries.  A  thermal  model  of  the  melt/rail  interface  is  used  to  determine  whether 
solidification  occurs  and  at  what  rate. 

An  additional  assumption  we  make  is  that  acceleration  is  insignificant.  In  the  experiments,  the 
sliders  were  accelerated  at  10  m/s  [8].  The  principal  effect  of  acceleration  is  to  produce  a  “lag” 
in  the  Couette  profile,  that  is,  there  is  a  time  required  for  the  velocity  distribution  to  adjust  to  the 
changing  speed.  The  lag  is  on  the  order  of  a  few  microseconds,  during  which  the  velocity 
changes  by  a  fraction  of  a  percent.  Acceleration  can  thus  be  neglected  because  the  velocity 
profile  is  lagging  by  only  a  small  fraction  of  the  velocity. 

Our  derivation  follows  closely  that  of  Stiffler  and  the  basic  equations  for  thin  films.  Following 
Stiffler,  variables  and  parameters  pertaining  to  the  liquid  film  are  not  subscripted  whereas  those 
pertaining  to  the  slider  and  rail  are. 
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MOMENTUM  AND  CONTINUITY  EQUATIONS 

The  momentum  equations  for  the  x  and  z  directions  are  (steady  state,  with  —  =  0 ) 

dy 

du  du  du  dp  d  du 

dx  dy  dz  dx  dy  dy  '  ' 

dw  dw  dv  dp  d  dw 

pu—+pv-—+pw—  +  -^=—(n—)  n\ 

dx  dy  dz  dz  dy  dy  ^ 

Assume  u=i^(y)  +  u2(x,y,z), with  M»|k2|,  where  Ul  is  the  Couette  flow  velocity  with  injection 

and  u2is  the  squeeze  film”  component  of  velocity.  Then,  applying  equation  (1)  to  u,  and  u2 

du 1  dii  dii  d  ,  dii 

PU'^+PV^  +  pW~A=^{,‘V  (3) 


du ,  du0  du ,  dp  d  du„ 

pu—1+pv-^+pw—l  +  ^=—(u^1):  (4) 

dx  ^  dy  H  dz  dx  dy*1  dy  {) 

Note  that  equation  (3)  describes  conditions  for  Couette  flow  with  —  =  0 .  Thus 

dx 

ui=uo@y=h 

Ui=0@y-0 

Also,  since  Uy  is  independent  of  x  and  z,  and  p,  is  assumed  constant,  equation  (3)  simplifies  to 

du,  d2u. 

<5> 

For  thin  films,  equations  (2)  and  (4)  can  be  simplified  by  neglecting  all  terms  on  the  left  hand 
. ,  dp  dp 

side  except  —  and  — .  Both  w  and  u2  are  assumed  to  be  independent  of  x  and  the  pressure 

gradients  due  to  momentum  convection  are  small  compared  to  the  viscous  forces.  Then  (2)  and 
(4)  become,  respectively 

dp  d2w 

dz  ~  P  df  (6) 


dp  _  d2u2 
dx'*  dy2 
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which  can  be  integrated  to  obtain 

1,2  ,  ,dp 

-yh,7*  (8) 

1,2  ,  .9/? 

~yh>Tz  (9) 

At  this  point,  it  is  convenient  to  introduce  a  new  variable  £,  the  fraction  of  the  melted  material 
that  solidifies.  (1-  Q  is  the  fraction  that  leaves  at  the  sides  of  the  slider.  Equations  8  and  9  show 
that  the  squeeze  film  flow  has  a  parabolic  profile  in  both  lateral  directions.  Applying  the 
continuity  equation  to  the  melt  film  at  distance  y  from  the  rail  (or  solidified  interface)  with  the 
overall  fraction  1  -  £  leaving  all  four  sides  with  the  parabolic  profiles  for  u2  and  w  gives 

v=v0(£ +(1-0(3 -2^hf)  (10) 

where  v0  =  .  The  fraction  £  is  determined  from  the  heat  transfer  analysis.  Then 

1  3d 

u  =  u1(y)  +  u2(  x,y,z)=ul(y)  +  —  (y2- yh)-f  (11) 

2(0.  dx 

where  is  obtained  from  numerical  integration  of  equation  (5).  The  solution  of  equation  (5)  is 
obtained  based  on  /^constant,  and  it  is  assumed  that  v  is  a  function  of  y  only. 


THE  FILM  REYNOLDS  EQUATION 

The  continuity  equation  for  the  film  is 

|(pu)+|(pv)  +  |(^)  =  0 

Substituting  u,v,  and  w  from  above  and  integrating  over  the  height  of  the  film  gives 


9  -p/i3  dp  f  d  -p/?3  9p  -9 
dx  12(0.  9x  9z  12 p.  dz  9x 


)=:~£[  fp^/+pv0a-o 


For  a  uniform  thickness  film, 


V2n  _  ~12mn(1-^ 
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This  is  the  same  relation  as  found  by  Stiffler,  [9]  except  for  the  l-C,  factor  which  accounts  for 
solidification.  The  film  thickness  is  given  by 


where  the  average  pressure  parameter,  8,  is  [10] 

<7=mean  normal  pressure  on  film  =  Fz/A 
Ls=smaller  ofL^  andLz 
192  00 

1—  ~5/3  S""5tanh(«^?/2)  (1 

"  P  rt=l,3,5 

P  =  Ratio  of  the  longer  side  of  the  slider  to  the  shorter  side  of  the  slider. 


DIMENSIONLESS  VARIABLES 

The  equations  are  rewritten  in  nondimensional  form  to  reduce  the  number  of  parameters  and 
obtain  solutions  in  terms  of  classical  nondimensional  quantities.  The  x  momentum  equation, 
equation  (5),  is  made  dimensionless  by 

v=V’  Y=y/h-  (17) 

^0  uo 

Because  Ux  is  only  a  function  of  Y,  the  result  is  written  as  an  ordinary  differential  equation, 

d2U,  dU, 

“7^L-Re^-Tf  =  0  (18) 


where  Rew  =  /?^.. 

A 


The  Y  velocity  can  be  rewritten  as 


V  =  -(^+(l-0(3r-2F3)) 

The  first  integration  is  performed  using  an  integration  factor  to  give 


-q  e+loRe»vdr 


where  Cu  is  a  constant  to  be  determined  by  the  boundary  conditions,  namely, 

c/,  =  o@y=o 
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Ul  =  l@Y=l  (23) 

Equation  (21)  was  integrated  numerically  using  a  fourth  order  Runge-Kutta  method.  The 
resulting  values  of  U^Y)  and  — 1  are  used  in  the  energy  equation. 


ENERGY  EQUATION 

The  boundary  layer  formulation  of  the  steady  state  energy  equation  for  the  liquid  film  with 
constant  specific  heat  and  thermal  conductivity  is  [  1 1  ] . 


dt  dt  dt  dp 

pcu—  +  pcv—  +  pew—  =  kVt  +  u—  +  utb 
ox  dy  dz  dx 


(24) 


where  the  dissipation  function  <p  is 


#-(1, 2+ 

dy  dy 


(25) 


Because  the  temperature  is  assumed  to  be  a  function  of  y  only  and  both  u  —  and  (— f  can  be 

dx  dy 

neglected,  the  energy  equation  simplifies  to 


dt  1  cPt  ,du* 

(26) 

We  make  the  energy  equation  dimensionless  using  equation  (17)  and 

6=--\ 

tm 

m 

(27) 

where  tm  is  the  melting  temperature.  The  result  is 

(28) 

where  u2  has  been  neglected  as  compared  to  and  where 

-aL 

ctm 

(29) 

Pr  =  — 
k 

(30) 
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BOUNDARY  CONDITIONS  FOR  THE  ENERGY  EQUATION 

A  total  of  four  boundary  conditions  must  be  matched  when  solving  the  liquid  film  energy 
equation.  These  are: 

1.  6(1)  -  0 .  The  liquid/  slider  interface  must  be  at  the  melting  point. 

2.  The  heat  flux  from  the  liquid  must  supply  the  melt  rate  energy  and  heat  conduction  to  the 
slider. 

3.  With  solidification  the  liquid  interface  must  be  at  the  melting  temperature,  0(O)=O. 
Without  solidification,  the  film  temperature  must  equal  that  of  the  rail  and  0(0)  must  be 
greater  than  zero. 

4.  With  solidification,  the  heat  flux  from  the  liquid  plus  the  solidification  energy  flux  must 
match  heat  conduction  to  the  resolidified  layer.  Without  solidification,  the  heat  flux  from 
the  film  must  match  the  flux  to  the  rail. 

Boundary  Condition  2 

Following  Stiffler  [9],  it  is  assumed  that  the  melting  rate  at  the  slider  has  reached  its  steady  state 
rate  where  the  flux  is  given  by  [12] 

qm=pv0[Hm+cs(tm-t0)]  (31) 

where  Hm  is  the  heat  of  fusion.  Thus,  the  dimensionless  heat  flux  boundary  condition  at  Y=1 
becomes 


[^fk,=-Re.Pr[^+(l-^)££] 

dY  ctm  t  c 


(32) 


The  boundary  condition  at  the  rail  side  of  the  liquid  film  is  more  complex.  The  two  cases, 
solidification  and  no  solidification,  are  considered  separately. 


No  Solidification  on  Rail  under  Slider 

The  temperature  and  heat  fluxes  at  the  rail/liquid  interface  must  match  and  the  interface 
temperature  must  be  above  the  melting  point.  The  heat  flux  to  the  rail  is  taken  to  be  the  average 
over  the  contact  time.  The  heat  flux  is  obtained  for  the  case  when  the  rail  surface  undergoes  a 
step  change  to  the  temperature  of  the  film  interface.  From  Carslaw  and  Jaeger  [13],  the  average 
heat  flux  into  a  semi-infinite  solid  from  t  =  0  to  t  =  Lx/u0  for  step  change  from  tr  to  tr(0)  is 


2kr(tr(0)-tJ 
\  uo 


(33) 


The  heat  flux  boundary  condition  at  the  rail/liquid  interface  is  then 
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dQ  _2k 
dY  r=0~  k 


where  f,(0)  =  [l  +  8(0)]tm. 


SOLIDIFICATION  ON  RAIL  UNDER  THE  SLIDER  v 

An  additional  interface  is  formed  when  solidification  occurs,  but  its  boundary  conditions  are 
included  in  the  solidification  solution.  The  approach  is  similar  to  the  case  for  no  solidification 
but  uses  the  solution  for  a  liquid  which  is  above  the  melting  temperature  that  comes  into  contact 
with  a  cold  solid  (see  Appendix  A). 


METHOD  OF  SOLUTION 

Once  one  has  assumed  a  value  for  vmJ  and  C,  one  can  calculate  the  film  thickness  for  the  given 
slider  dimensions,  pressure,  and  velocity.  One  can  also  integrate  the  x  momentum  equation  to 

obtain  values  of  Uj  and  .  Because  the  two  interface  conditions  at  the  slider,  conditions  (1) 

and  (2),  are  also  known,  one  can  integrate  the  second  order  energy  equation  starting  at  Y=1  and 
proceeding  to  Y=0.  Note  that  both  initial  value  and  slope  are  available  at  Y=l.  The  result  of  the 

integration  will  be  for  0(0)  and  ~~  .  The  iteration  is  completed  first  for  no  solidification. 

For  no  solidification,  one  keeps  £  =  0  and  changes  vms  until  the  heat  flux  from  the  film  matches 
the  heat  flux  to  the  rail  for  the  value  of  0(0).  If  the  value  of  0(0)  is  greater  than  0  when  the  fluxes 
match,  there  is  no  solidification  and  the  solution  is  finished.  Otherwise,  one  begins  the  iteration 
for  solidification. 

For  solidification  on  the  rail,  one  must  first  find  a  non-zero  value  of  £  needed  to  satisfy  the 
boundary  condition  (3),  #(0)  =  0.  If  #(0)is  not  zero,  then  the  assumed  value  of  C,  is  adjusted 
until  it  is  zero.  When  this  has  been  accomplished,  only  boundary  condition  (4)  remains  to  be 
satisfied.  A  value  of  tv  is  assumed  and  X  is  obtained  from  the  solidification  solution  along  with 
the  values  of  X  and  of  the  heat  flux  at  the  interface.  The  value  of  tv  is  adjusted  until  the  flux  at 
the  solidified  surface/film  interface  matches  the  flux  from  the  film  solution.  When  the  heat 
fluxes  at  the  solidified  boundary  match,  it  is  then  determined  if  the  value  of  X  is  too  high  or  too 
low  to  match  the  assumed  solidification  fraction,  £.  If  not,  vm  is  adjusted  accordingly,  and  the 

procedure  is  repeated  until  the  solidification  rates  also  match. 


RESULTS 

The  model  is  run  with  the  objective  of  reproducing  the  experimental  results  described  in  the 
introduction  of  this  paper.  Three  slider  pressures  are  considered  10.9,  17.0,  and  21.8  ksi,  which 
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effectively  bracket  the  data  at  pressures  characteristic  of  railgun  armatures.  For  each  value  of 
pressure,  solutions  are  found  for  slider  speeds  ranging  from  800  m/s  to  1800  m/s  using  the 
parameters  and  properties  shown  in  Table  1.  We  used  the  incipient  melting  temperature,  805  K, 
as  the  melting  temperature. 


Symbol 

Variable  Name 

Table  1. 

Physical  Quantity 

Value 

Units 

p 

RHOL 

Density  of  liquid 

2485 

kg/m3 

K 

CONDKL 

Thermal  conductivity  of  liquid 

83.7 

W/m/K 

ii 

MU 

Absolute  viscosity  of  liquid 

0.00450 

Ns/m2 

Pr 

RHOR 

Density  of  rail 

8950 

kg/m3 

K 

CONDKR 

Thermal  conductivity  of  rail 

374 

W/m/K 

Ps 

RHOS 

Density  of  slider 

2800 

kg/m3 

K, 

CONDKS 

Thermal  conductivity  of  slider 

170 

W/m/K 

C 

CPL 

Specific  heat  of  liquid 

1084 

J/kg/K 

cr 

CPR 

Specific  heat  of  rail 

400 

J/kg/K 

c. 

CPS 

Specific  heat  of  slider 

1200 

J/kg/K 

HM 

Latent  heat  of  melting  of  slider 

3.78xl05 

J/kg 

4 

LZ 

Length  of  slider 

0.0155 

m 

h 

LX 

Width  of  slider 

0.0014 

m 

TINSM 

Incipient  melting  temperature  of  slider 

805 

K 

TLIQ 

Liquidus  temperature  of  slider 

908 

K 

Jm _ 

TM 

Melting  temperature  of  slider 

805 

K 

MELT  VELOCITY,  ORIGINAL  MODEL 

A  direct  comparison  between  model  and  experiment  is  obtained  by  comparing  Figure  3,  which 
plots  calculated  melt  velocity  vms  vs.  slider  speed,  with  experimental  curves  in  Figure  1. 
Although  the  model  predicts  melt  velocities  that  are  centered  within  the  range  of  experimental 
values  (vms=  1.1  m/s),  the  model  fails  to  reproduce  the  strong  dependence  on  velocity  and 
pressure  observed  in  the  experiments. 
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Melt  Velocity  vs.  Slider  Velocity 


Skier  Velocity,  nVa 


Fig.  3.  Initial  results  for  three  values  of  pressure.  The  relative  insensitivity  of  the  model  to  velocity  and  pressure 
shows  poor  agreement  with  experiment. 

The  relative  insensitivity  of  the  model  to  velocity  and  pressure  can  be  attributed  to  the 
solidification  boundary  condition  (4)  coupled  with  the  fact  that  residence  time  is  inversely 
proportional  to  the  speed  of  the  slider.  Although  there  is  an  increased  average  heat  transfer  rate 
(transient  conduction  into  the  rail),  there  is  less  solidification  under  the  slider  as  speed  increases. 
Reduced  solidification  at  high  speeds  means  that  the  film  has  to  get  thicker,  which  prevents 
viscous  dissipation  from  increasing  appreciably.  This  strong  self-limiting  effect  is  a  feature  of 
the  laminar  model  with  high  solidification  rates.  The  fraction  of  melt  solidified  under  the  slider 
ranged  from  81.5  to  99.7  percent  over  the  range  of  conditions  in  Figure  3. 

The  small  temperature  rise  in  the  film  also  indicated  a  problem  with  the  model  when  we 
reconsidered  the  large  range  between  the  incipient  melting  temperature  and  the  liquidus 
temperature.  The  peak  temperature  in  the  film  of  817  K,  which  is  12  K  above  the  incipient 
melting  temperature  and  92  K  below  the  liquidus  temperature,  is  inconsistent  with  the 
assumption  that  the  full  latent  energy  has  been  absorbed  by  the  molten  aluminum. 


REVISED  MODEL 

This  realization  motivated  us  to  include  a  more  realistic  approximation  to  the  complex  phase 
change  behavior  of  7075  aluminum  alloy.  We  also  reasoned  that  lowering  the  energy  needed  to 
release  aluminum  into  the  film  ought  to  increase  the  melt  speed  vms.  Additionally,  we  reasoned 
that  the  viscosity  near  the  incipient  melting  temperature  should  be  much  higher  than  that  of 
liquid  aluminum.  Appendix  B  describes  our  effort  to  incorporate  a  description  of  distributed 
latent  heat  into  our  model.  As  a  brief  summary  of  the  revised  model,  we  kept  approximately  10 
percent  of  the  latent  heat  at  a  temperature  of  815  K,  10  K  above  the  incipient  melting  point.  The 
remainder  of  the  latent  heat  was  assumed  uniformly  distributed  over  the  range  from  815  to  908 
K,  the  liquidus  temperature. 
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We  also  attempted  to  develop  an  appropriate  model  of  the  viscosity  of  the  partially  liquid  state  of 
7075  aluminum  between  805  and  908  K.  The  variable  viscosity  model  is  also  described  in 
Appendix  B.  The  model  was  ultimately  discarded  because  the  effect  of  decreasing  viscosity  with 
increasing  temperature  predicted  a  smaller  increase  in  melt  velocity  with  slider  speed  than 
resulted  from  a  constant  viscosity. 


The  melt  velocity  results  for  the  revised  model  are  shown  in  Figure  4.  As  expected,  the  rate  of 
melting  was  increased,  but  the  rates  were  still  not  very  sensitive  to  pressure  or  slider  velocity. 


Melt  Velocity  vs.  Slider  Velocity 


Fig.  4.  Results  of  model  revised  to  include  a  more  realistic  approximation  to  the  phase  change  behavior  of  7075 
aluminum  alloy  and  a  model  of  temperature-dependent  model  of  viscosity.  As  expected,  the  rate  of  melting 
increased  but  the  rates  are  not  very  sensitive  to  pressure  or  slider  velocity. 


MODEL  INADEQUACIES 

The  poor  agreement  indicates  that  one  or  more  of  the  assumptions  made  in  the  model  is 
incorrect.  We  examined  several  of  the  assumptions  that  appeared  to  be  most  questionable.  These 
include:  (1)  perfect  thermal  contact  between  film  and  rail,  (2)  laminar  flow  of  the  film,  and  (3) 
simple  melting  of  the  aluminum  alloy. 

We  first  explored  the  effect  of  a  thermal  interface  resistance  between  the  film  and  the  rail.  Early 
experimental  work  on  liquid  metal  heat  transfer  was  unable  to  confirm  the  high  heat  transfer 
rates  predicted  at  low  Prandtl  numbers.  It  turned  out  that  interface  resistance  due  to  very  thin 
oxide  layers  and  non-wetting  greatly  reduced  the  heat  transfer  from  the  theoretical  rates.  The 
same  effect  might  be  present  in  the  railgun  environment.  It  could  be  caused  by  similar 
phenomena  or  by  gases  entrained  as  the  slider  advances. 

Because  our  formulation  could  not  readily  be  modified  to  include  this  in  the  boundary  condition, 
we  approximated  the  effect  by  reducing  the  thermal  conductivity  of  the  rail.  A  reduction  of 
thermal  conductivity  of  the  rail  by  a  factor  of  ten  results  in  a  reduction  of  heat  transfer  to  the  rail 
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by  about  a  factor  of  three.  Results  still  exhibited  the  same  trends  as  before  for  the  nominal  liquid 
viscosity. 

Another  issue  that  arose  was  the  validity  of  the  laminar  flow  assumption.  Early  calculations 
using  Stiffler’s  model  without  solidification  indicated  that  the  film  could  be  expected  to  be  in  a 
turbulent  flow  regime.  In  most  of  our  program  data  runs,  the  Reynolds  number  for  the  film 
reached  the  turbulent  transition  at  the  higher  slider  velocities.  We  decided  to  approximate  the 
effects  of  turbulence  by  increasing  the  eddy  diffusivity  of  heat  and  momentum.  We  applied  a 
factor  to  increase  both  the  viscosity  and  thermal  conductivity  by  a  uniform  amount  across  the 
width  of  the  channel.  When  we  combined  this  with  the  reduced  heat  transfer  to  the  rail,  we 
obtained  results  that  closely  followed  the  major  trends  of  the  test  data. 

Examples  of  the  results  are  shown  in  Figures  5  and  6.  In  Figure  5,  the  rail  thermal  conductivity 
is  reduced  by  a  factor  of  ten  and  the  viscosity  and  thermal  conductivity  of  the  melt  film  are  each 
increased  by  a  factor  of  three.  The  curve  for  a  =  10.9  ksi  agrees  very  well  with  the  test  data  in 
Figure  1  between  1300  and  1800  m/s.  Figure  6  is  for  the  conditions  with  rail  thermal 
conductivity  reduced  by  a  factor  of  15,  and  both  viscosity  and  thermal  conductivity  of  the  melt 
increased  by  a  factor  of  four.  The  curve  in  Figure  6  for  a  =  21.8  ksi  is  in  very  good  agreement 
with  the  test  data  from  1000  to  1800  m/s. 


Melt  Velocity  vs.  Slider  Velocity 


Sfider  Velocity,  m/s 


Fig.  5.  Results  of  model  revised  to  (1)  suppress  heat  transfer  to  the  rails  by  reducing  the  thermal  conductivity  of  the 
rails  by  a  factor  of  10  and,  (2)  approximate  turbulence  by  increasing  viscosity  and  conductivity  by  a  factor  of  three. 
The  change  in  slope  between  1100  and  1300  m/s  corresponds  to  the  end  of  melt  solidification  under  the  slider. 
Above  1300  m/s,  the  trends  agree  well  with  the  data. 
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Melt  Velocity  vs.  Slider  Velocity 


Fig.  6.  Results  of  the  model  revised  to  (1)  suppress  heat  transfer  to  the  rails  by  reducing  the  thermal  conductivity  of 
the  rails  by  a  factor  of  15  and,  (2)  approximate  turbulence  by  increasing  viscosity  and  conductivity  by  a  factor  of 
four.  The  changes  in  slope  at  900  to  1000  m/s  correspond  to  the  end  of  melt  solidification  under  slider.  Good 
agreement  with  experiments  above  the  changes  in  slope  suggests  that  the  film  does  not  solidify  under  the  armature. 

A  change  in  slope  is  seen  in  Figure  5  at  slider  velocities  between  1 100  and  1300  m/s.  Above  this 
change  in  slope,  the  trend  agrees  well  with  the  test  data.  In  Figure  6,  the  slope  changes  occur 
between  900  and  1000  m/s.  Again  the  trends  agree  well  above  the  change.  The  change  in  slope 
corresponds  to  the  end  of  melt  solidification  under  the  slider,  suggesting  that  the  film  does  not 
solidify  under  the  armature  at  all. 


DISCUSSION 

There  are  several  observations  that  can  be  made. 

1.  The  model  results  do  not  apply  for  slider  velocities  below  approximately  1000  m/s.  This  is 
believed  to  be  due  to  the  transient  nature  of  the  railgun-driven  wear  experiments.  In  a 
separate  analysis  of  solid-solid  sliding,  we  found  that  the  melting  temperature  is  reached  at 
about  the  time  the  velocity  reaches  700  m/s  for  an  acceleration  of  106  m/s2.  Thus,  we  do  not 
expect  the  assumption  of  quasi-steady-state  melt  to  be  satisfied  until  a  velocity  of  1000  m/s  is 
reached. 

2.  Reduced  heat  transfer  to  the  rail  combined  with  higher  viscosity  and  thermal  conductivity  of 
the  melt  film  provides  a  better  fit  to  the  data.  A  reduction  in  rail  heat  transfer  due  to  interface 
thermal  resistance  of  a  factor  of  three  to  four  appears  to  be  needed.  Effective  viscosity  and 
thermal  conductivity  higher  than  their  nominal  liquid  values  by  a  factor  of  three  to  four  also 
appear  to  be  needed.  Viscosity  could  be  explained  either  by  turbulence  or  by  incomplete 
melting  because  the  maximum  liquid  film  temperature  was  below  the  liquidus  temperature 
for  all  conditions  in  Figures  5  and  6. 
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3.  The  results  are  not  highly  sensitive  to  the  location  of  the  pseudo  melting  point,  so  long  as  it  is 
within  20  percent  of  the  incipient  melting  temperature. 

The  results  give  encouragement  that  a  valid  model  can  be  developed.  The  model  clearly  indicates 
that  film  solidification  under  the  armature  is  not  an  important  process.  Future  work  should  focus 
on  modifying  the  analysis  to  include: 

•  turbulent  correlations  of  eddy  viscosity  and  eddy  diffusivity, 

•  additional  information  on  viscosity  and  the  distribution  of  heat  of  transformation  in  the 
range  between  the  incipient  melting  and  liquidus  temperatures, 

•  the  effect  of  thermal  interface  resistance  between  the  film  and  the  rail. 

The  complexities  of  resolidification  can  be  avoided  in  future  work. 
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APPENDIX  A.  THERMAL  BOUNDARY  CONDITIONS  FOR  SOLIDIFICATION  ON 

THE  RAIL 


The  solidification  boundary  conditions  are  satisfied  by  applying  the  results  of  a  classical 
solidification  problem:  the  region  x  <  0  initially  liquid  and  the  region  x<0  solid  [14]. 

The  region  x  <  0  has  thermal  constants  pr,cr,kr,Kr  and  is  initially  at  temperature  tfo .  The 
region  x  <  0  is  initially  liquid  with  thermal  constants  p,c,k,K  and  is  initially  at  temperature  tv. 
The  thermal  constants  ps,cs,ks,Ks  of  the  solidified  liquid  in  the  region  0  <  x  <X  may  be 
different  than  the  liquid. 

The  temperatures  at  time  x  of  the  solid  (rail),  the  solidified  layer  and  the  liquid,  tr,  ts  and  t, 
respectively,  are  given  by: 

tr  =  AA(l  +  erf(— f=))  +  tr \x<0 

2yjJCrT 

(A-l) 

t,  =BB  +  CCerf  (— *  )  +  t',0<x<X 

(A-2) 

t  =  tv-  DDerfc{ — j=)‘,  x>X 

2-slfcr 

(A-3) 

The  position  of  the  interface  at  time  t  is 

X  =  2/1  yjxsT 

(A-4) 

where 

exp(-i)  BZexp(-MSL) 

[BSR  +  erf(X)  BSLerfc(AASL) 

(A-5) 

AA=BB=BSR  (tm-tJ/(BSR  +  erf(X)) 

(A-6) 

CC  =  AA/ 

^  /BSR 

(A-7) 

DD  =(tv  -  tm)  /  erfc(A  •  ASL) 

(A-8) 

ASI.=  -K' 

(A-9) 

BZ =(tv-tm)/(t.-tj 

(A- 10) 
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(A-ll) 


\ 


The  process  of  solidification  of  the  melt  film  is  considered  equivalent  to  the  above  solidification 
process  when  the  average  value  of  the  gradient  at  the  melt  film  matches  the  above  value  and  the 
thickness  of  the  solidified  layer  also  matches  the  thickness  assumed  for  the  melt  film 
calculations.  The  initial  liquid  temperature,  tv,  is  used  as  a  parameter  to  match  the  gradient  at  the 
interface  in  the  two  solutions. 
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APPENDIX  B.  REVISED  LATENT  HEAT  AND  VISCOSITY  MODEL 

An  alloy  that  exhibits  an  incipient  melting  and  a  liquidus  temperature  undergoes  a  phase  change 
in  which  the  energy  of  transformation  is  distributed  over  the  range  between  the  two 
temperatures.  The  distributed  transformation  from  solid  to  liquid  is  also  accompanied  by  a 
change  from  the  infinite  viscosity  at  incipient  melting  to  the  viscosity  of  the  liquid  at  the  liquidus 
temperature. 

Because  information  on  the  distribution  of  the  energy  is  not  available,  we  assumed  that  it  occurs 
uniformly.  The  thermal  behavior  during  this  melting  process  can  be  modeled  by  adding  the 
latent  heat  of  melting  per  degree  temperature  change  to  the  normal  specific  heat.  The  extra 
specific  heat  applies  to  both  heating  and  cooling  within  the  range  of  transformation.  We 
arbitrarily  chose  the  pseudo  melting  point  to  be  815  K,  10  K  above  the  incipient  melting 
temperature.  The  fraction  (815-805)/(908-805)  =  0.0971  of  the  latent  heat  was  assumed  to  occur 
at  this  pseudo  melting  point.  The  remaining  fraction,  0.9029,  was  distributed  uniformly  as  an 
increased  specific  heat  between  815  and  908  K.  For  a  latent  heat  of  3.78  X  10s  J/kg,  the  increase 
in  liquid  specific  heat  is  3670.  J/kg/K  over  this  temperature  range. 

The  viscous  behavior  of  the  material  in  this  range  is  more  complex.  To  greatly  simplify  the 
problem,  it  was  assumed  that  the  fraction  of  the  energy  of  transformation  corresponds  to  the  local 
fraction  of  liquid,  then  the  viscous  behavior  can  be  related  to  a  two-phase  mixture.  The  simplest 
model  is  that  of  micro  layers  of  liquid  and  solid,  the  thickness  of  the  liquid  being  $  of  the  total, 
where  <j)  is  the  liquid  fraction.  The  local  viscosity  for  this  configuration  is  given  by  jxm  = 

In  order  to  retain  the  basic  formulation,  we  assumed  that  for  §  <  (j)0  the  viscosity  remains  infinite 
and  that  the  fraction  <J>0  of  the  latent  heat  is  concentrated  at  the  corresponding  pseudo  melting 
temperature,  tm'=  W  +(tL-tim)(J)o  where  tim  is  the  incipient  melting  temperature  and  tL  is  the 
liquidus  temperature.  The  fraction  l-4>o  of  the  latent  heat  is  uniformly  distributed  between  tm' 
and  tL.  Note  that  for  tm'  =  815  K,  <J>o  =  0.0971.  The  modified  model  applies  as  long  as  the 
maximum  temperature  in  the  melt  film  is  below  t^ 

This  model  of  the  phase  change  process  was  included  in  the  squeeze  film  analysis  and  in  the 
thermal  energy  equation  by  finding  effective  viscosity  for  a  film  of  thickness  h  with  a  parabolic 
temperature  profile  using  the  local  viscosity  model  described  above.  The  effective  viscosity  for 
the  squeeze  film  flow  is  obtained  from  the  solution  of  the  x  direction  momentum  equation 
neglecting  the  effect  of  injection  and  suction.  The  value  of  4>  can  be  expressed  in  terms  of  y  for  a 
parabolic  temperature  profile,  0  =  40cY(l-Y)  by 

where  b  =  40Ctm'/(tL- tim).  Then  the  momentum  equation 


is  integrated  for  constant 


— -  using  (i  =  (j.L/({)  with  u  =  0  at  y  =  0  and  u  =  0  at  y  =  h. 
oz 


Comparing 


the  result  with  the  solution  to  the  same  equation  for  p  =  Pl  gives 
M-s  =  |iL/(<t>0+b/10) 

where  ps  is  the  effective  viscosity  for  calculating  the  value  of  h  for  squeeze  flow  using  the 
analysis  for  uniform  viscosity. 


A  similar  process  is  repeated  to  find  the  solution  for  Couette  flow  with  the  same  value  of  <j>. 
After  one  obtains  the  velocity  gradient,  one  integrates  to  find  the  total  dissipation.  Comparing 
the  result  with  the  same  process  for  |1  =  pL  gives 
Pq  =  |i/(<j>o  +  b/6) 

where  Pq  is  the  effective  viscosity  used  to  calculate  the  viscous  heating  for  Couette  flow  using  a 
constant  viscosity  analysis. 
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